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ABSTRACT 

We study cluster mass dark matter haloes, their progenitors and surroundings in an coupled 
Dark Matter-Dark Energy model and compare it to quintessence and ACDM models with adi¬ 
abatic zoom simulations. When comparing cosmologies with different expansions histories, 
growth functions & power spectra, care must be taken to identify unambiguous signatures of 
alternative cosmologies. Shared cosmological parameters, such as ag, need not be the same 
for optimal fits to observational data. We choose to set our parameters to ACDM 2 : = 0 values. 
We find that in coupled models, where DM decays into DE, haloes appear remarkably similar 
to ACDM haloes despite DM experiencing an additional frictional force. Density profiles are 
not systematically different and the subhalo populations have similar mass, spin, and spatial 
distributions, although (sub)haloes are less concentrated on average in coupled cosmologies. 
However, given the scatter in related observables (14iax, ^y^ax)’ this difference is unlikely to 
distinguish between coupled and uncoupled DM. Observations of satellites of MW and M31 
indicate a significant subpopulation reside in a plane. Coupled models do produce planar ar¬ 
rangements of satellites of higher statistical significance than ACDM models, however, in all 
models these planes are dynamically unstable. In general, the nonlinear dynamics within and 
near large haloes masks the effects of a coupled dark sector. The sole environmental signa¬ 
ture we find is that small haloes residing in the outskirts are more deficient in baryons than 
their ACDM counterparts. The lack of a pronounced signal for a coupled dark sector strongly 
suggests that such a phenomena would be effectively hidden from view. 

Key words: (cosmology:) dark matter, (cosmology:) dark energy, galaxies:clusters:general, 
methods: numerical 


1 INTRODUCTION 

The ACDM concordance model of cosmology describes a spatially 
flat universe with an energy budget dominated by a invisible Dark 
Sector composed of two major components: Dark Matter (DM) that 
governs the small scale clustering of luminous matter; and Dark En¬ 
ergy (DE), which drives the late time accelerated expansion. This 
model is based on multiple lines of evidence, from the Cosmic Mi¬ 
crowave Background (CMB) anisotropies (e.g. Bennett et al. 2013; 
Planck Collaboration et al. 2013, 2015), Baryonic Acoustic Oscil¬ 
lations (BAO) (Beutler et al. 2011; Blake et al. 2011; Anderson 
et al. 2014, e.g), Large-Scale Structure (LSS) (e.g. Abazajian et al. 
2009; Beutler et al. 2012), weak lensing (e.g. Kilbinger et al. 2013; 
Heymans et al. 2013), cluster abundances (e.g Vikhlinin et al. 2009; 
Rozo et al. 2010), galaxy clustering (e.g. Tegmark et al. 2004; Reid 
et al. 2010), to the luminosity-distance relation from Type la (e.g. 
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Kowalski et al. 2008; Conley et al. 2011; Suzuki et al. 2012). De¬ 
spite all this evidence, the nature of the Dark Sector remains one 
of the most fundamental mysteries in cosmology. Observational 
evidence strongly favours nonbaryonic elementary particle DM, 
though the exact properties of the DM particle(s) are poorly con¬ 
strained. Most studies have focused on so-called Cold Dark Mat¬ 
ter (CDM) (see Frenk & White 2012, for a review), which has 
well-motivated candidates from particle physics, e.g. the lightest 
supersymmetric particle or neutralino (e.g. Ellis et al. 1984; or for 
a summary of several candidates see Bertone et al. 2005; Petraki & 
Volkas 2013). Similarly, the measured expansion rate at late times 
implies that the DE is well characterised by a cosmological con¬ 
stant, A, with an measured equation-of-state (eos) w = p/p ~ — 1, 
although the fact that the Universe is not exactly homogeneous may 
introduce bias in these observations (e.g. Clarkson et al. 2012). 

In spite of the remarkable agreement of ACDM with obser¬ 
vations of large-scale structure, the model is in tension with ob¬ 
servations on galaxy scales, the most well known of which is the 
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so-called “missing satellite problem”: CDM models predict many 
more satellite galaxies than observed around galaxies such as our 
own (e.g. Klypin et al. 1999; Moore et al. 1999). The excess sub¬ 
haloes that do not host galaxies may indicate that feedback pro¬ 
cesses, such as surpernovae, efficiently remove gas from low mass 
subhaloes leaving them almost completely dark (e.g. Bullock et al. 
2000; Benson et al. 2002; Nickerson et al. 2011, 2012), or may 
indicate the need for modifications to CDM, such as Warm Dark 
Matter (e.g. Lovell et al. 2014; Schneider et al. 2014; Power 2013; 
Elahi et al. 2014). Perhaps the most intriguing observational co¬ 
nundrum is the as yet unexplained alignments of satellite galaxies 
surrounding the Milky Way and its galactic neighbour, Andromeda, 
the so-called Vast Polar Structure (VPOS Pawlowski et al. 2012) 
and Plane of Satellites (Ibata et al. 2013; Conn et al. 2013) respec¬ 
tively. This alignment is not unique to the local group; alignments 
have observed in SDSS data (e.g. Yang et al. 2006; Li et al. 2013), 
and observations show satellite galaxies have velocities that are 
significantly more correlated than theory would predict (e.g., Ibata 
et al. 2014; Gillet et al. 2015; for counterarguments see Sawala et al. 
2014; Libeskind et al. 2015). 

There are also some theoretical shortcomings with ACDM, 
such as the fine-tuning and coincidence problems, which remain 
poorly explained from a purely theoretical point-of-view (see for 
instance Ferreira & Joyce 1997; Brax & Martin 2000; Huey & Wan- 
delt 2006). The former problem refers to the fact that, if we assume 
that the DE is a cosmological constant arising from the zero-point 
energy of a fundamental quantum field, then its density requires 
an unnatural fine-tuning of several tens of orders of magnitude to 
be compatible with observed cosmological constraints. The latter 
problem arises from the difficulty in satisfactorily explaining the 
observation that matter and DE energy densities today have com¬ 
parable values given that the energy density of these two fluids have 
completely different dependencies on cosmic time and only now is 
DE beginning to dominate. If DE dominated at early times, cosmic 
structure formation would be strongly suppressed. 

This coincidence problem has lead to the proposal of several 
alternatives, such as modifications to general relativity (e.g. Hu & 
Sawicki 2007; Starobinsky 1980; and de Felice & Tsujikawa 2010 
for a review) or dynamical scalar fields (e.g. Ratra & Peebles 1988; 
Wetterich 1988; Armendariz-Picon et al. 2001; and Tsujikawa 2013 
for a review). Dynamical scalar fields, or quintessence models, are 
an attractive alternative since, in principle, they can alleviate some 
of the fine-tuning problem. An interesting subset of these models 
are those where the scalar field couples to the matter sector, thereby 
removing the coincidence problem (e.g. Wetterich 1995; Amendola 
2000). Typically it is assumed that the dark energy couples only to 
dark matter, be it cold or warm, (for lack of evidence indicating 
normal matter interacting with a hidden sector). The coincidence 
problem is alleviated by having the dark matter decay to the scalar 
field producing the late time accelerated expansion. 

In N-Body realisations of these models, this interaction gives 
rise to two novel effects: the DM N-Body particle masses vary with 
time; and the effective gravity constant governing two-body inter¬ 
actions between DM N-Body particles is no longer the same as that 
governing baryon-DM or baryon-baryon interactions. Studies have 
examined cosmological structure formation using N-Body simu¬ 
lations for both uncoupled (e.g. Klypin et al. 2003; Dolag et al. 
2004; De Boni et al. 2011) and coupled models (e.g. Maccio et al. 
2004; Baldi 2012; Baldi et al. 2010; Li & Barrow 201 1). Li & Bar- 
row (201 1) showed that there are significant differences in the mat¬ 
ter power spectrum as a result of the different growth functions. 
Sutter et al. (2015) found void density profiles and the number of 


very large voids are affected by coupled cosmologies, although this 
study compared ACDM to a very strongly coupled model, one that 
is ruled out by observations (Pettorino et al. 2012) and the differ¬ 
ences are small. Pace et al. (2015) & Giocoli et al. (2015) found 
differences in the weak lensing signature of coupled models. How¬ 
ever, Carlesi et al. (2014a,b) found negligible differences in the cos¬ 
mic web and halo mass function between ACDM uncoupled and 
coupled models. However, these authors did find weak differences 
in the concentration and spin parameter of small (moderately re¬ 
solved) field haloes. 

The hint that (un)coupled quintessence models statistically 
differ from ACDM at not only large scales but for small haloes 
is exciting given the tensions that exist between observations and 
ACDM predictions. Our aim is to examine the distribution of dark 
matter haloes in hydrodynamical zoom simulations, identify the 
observational fingerprints of a coupled dark sector and explore 
whether these models can reduce the observational tensions that 
exist with the current concordance ACDM model. This paper is or¬ 
ganised as follows: we describe the numerical methods in §2 and 
discuss some issues concerning comparing different cosmological 
models. Our findings are presented in sections 3-5. These results 
are discussed in §6. 

2 METHODS 
2.1 The model 

As the theoretical and numerical basis for coupled quintessence 
models has been thoroughly covered before we only briefiy de¬ 
scribe the cosmological model here (for theoretical and observa¬ 
tional discussions see Amendola 2000; Tsujikawa 2013; for discus¬ 
sions of numerical approaches see for instance Maccio et al. 2004; 
Baldi 2012; Li & Barrow 2011; Carlesi et al. 2014a). Dark energy 
in these models arises from the evolution of a scalar field, 0, whose 
Lagrangian is generically written as: 

L = y* + y(0) + m(0)'0m'0m] , (1) 

that is a kinetic term, a potential term and a term characterising the 
coupling of the scalar field to the matter sector. There are several 
models for the potential, V (0) that give rise to late time accelerated 
expansion, mimicking the affect of a cosmological constant A. Here 
we focus on the so-called Ratra-Peebles potential (Ratra & Peebles 
1988): 

v{4>) = ( 2 ) 

where 0 is in units of the Planck mass. 

There is a great deal of freedom in choosing the form of m(0), 
the interaction term. A popular, simple model that can alleviate the 
coincidence problem is an exponential coupling of the scalar field 
to dark matter only 

m{(j)) = rrio exp [-/3(0)0] . (3) 

For simplicity, we use a constant coupling term, /3(0) = fio. The 
first consequence of this coupling is the decay of dark matter parti¬ 
cles into the scalar field. A second consequence is that dark matter 
and baryonic matter are governed by different dynamics. Whereas 
baryons follow standard Newtonian dynamics, dark matter experi¬ 
ences an additional frictional force which appears as an effective 
gravitational force, Geff = G(1 + /3o). The two matter fluids de¬ 
velop an offset in the amplitude of their density perturbations, a 
’gravitational bias’ that can significantly impact the baryon frac¬ 
tion of galactic- or cluster-sized haloes (e.g. Baldi et al. 2010). 
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2.1.1 Initial conditions and comparing cosmologies 

Before discussing the simulations themselves it is worthwhile con¬ 
sidering the initial conditions and how cosmological parameters 
such as the matter density are set. Different cosmological obser¬ 
vations tightly constrain different sets of cosmological parameters, 
with each observation having different correlations between vari¬ 
ous parameters. CMB measurements, for instance, primarily con¬ 
strain quantities such as the amplitude of the matter power spectrum 
at a moment in time during the matter dominated era when the last 
scattering surface is produced ^ . Supernovae observations measure 
the expansion rate over cosmic time at (moderately) low redshifts. 
Moreover, observational data is only meaningful in the context of a 
model. As a consequence, some models and their associated param¬ 
eters are not necessarily tightly constrained if one only considers a 
single observational data set. 

An alternative model may share cosmological parameters with 
ACDM, such as a8{z = 0), the normalisation of the mass variance 
quantifying the nonlinearity of regions enclosing a fixed amount 
of mass. However, analysing the data with a particular model will 
lead to different optimal values for these shared parameters, which 
may not even agree to within some confidence level. This gives 
some flexibility in setting the cosmological parameters when run¬ 
ning simulations of alternative models. One could choose to set the 
values of the shared subset based on estimates derived from the 
CMB interpreted through the ACDM lens. A consequence of this 
choice is that the parameters could have far different z — ^ values 
in the alternative model when compared to ACDM. An example 
of the above are the coupled models studied by Baldi (2012). This 
author fixed the power spectrum amplitude at zcmb, ^s, across all 
the models but Ob & Om are set to be equal at z = 0. The re¬ 
sult is that as differs significantly between models, from 0.8 — 0.9. 
In contrast, the matter density parameters start at different values 
in order to agree at z = 0 with a ACDM interpretation of CMB 
observations. 

Here we have taken a slightly different approach, used in Car- 
lesi et al. (2014a), and set both as{z — 0) and the matter density 
parameters based on CMB observations by Planck interpreted with 
a ACDM model (Planck Collaboration et al. 2013). As a conse¬ 
quence, the amplitude of the density perturbations matches at = 0 
but starts from a different point at ^cmb and zi, the redshift were 
we start our simulations. Neither approach is more correct than the 
other. However, it is important to note that differences will arise 
between cosmological models due to “simple” differences in cos¬ 
mological parameters like as and Vtrn that may hide differences 
arising from alternative physics or differences in the growth of den¬ 
sity perturbations ^ We note that although methods exist that map 
the results from a simulated cosmology to another, producing the 
same halo mass function, the resulting mapped haloes have sys¬ 
tematically biased internal properties, i.e., the nonlinear evolution 
is not correctly accounted for (e.g. Angulo & White 2010; Mead 
& Peacock 2014). Thus we do not attempt to estimate the change 
in a given property is a result of different cosmological parameters, 
rather we attempt to identify properties that appear affected by the 
inclusion of extra dark sector physics. 


^ This is not strictly true as distortions along the line-of-sight, from weak 
leasing for example, affects the temperature power spectrum of the CMB 
^ See Pace et al. (2015), who show the change in weak leasing signals from 
a ACDM model with the same as (z = 0) as an alternative model accounts 
for a significant portion of the signal in the alternative model. 


Table 1. Comparison of several key cosmological parameters relative 
ACDM at 2 ; = 0.5. Here ttidm is dark matter mass, D+ is the growth 
function, H is the Hubble expansion rate, and Geff is the effective gravita¬ 
tional potential in DM-DM interactions. 


Cosmology 


^^cdm 

D+ 

H 

Geff 

ACDM 

(-) ACDM 

“CDM 

7-) ACDM 

^ACDM 

^ACDM 

Ggff 

0CDM 

1 

0.979 

0.992 

1.0129 

1 

<f>(0 = 0.05)CDM 

1.0027 

0.993 

1.009 

1.0024 

1.0025 


2.2 Simulations 

We use 3 cosmologies in this study: a reference ACDM; an uncou¬ 
pled quintessence model ((/)CDM); and one coupled model (/lo = 
0.05). As mentioned previously all cosmologies have parameters 
consistent with z = 0 ACDM Planck data (h, Qm, ^b, os) = 
(0.67,0.3175,0.049,0.83) (Planck Collaboration et al. 2013, 
2015). Our coupled model is chosen to test the boundaries of al¬ 
lowed coupling (see Pettorino et al. 2012) in order to maximise 
any observational differences that may exist between this cosmol¬ 
ogy and the standard ACDM model. The power spectrum and 
evolution of these quantities along with the growth factor / = 
din D(a)/d In a are calculated using first-order Newtonian pertur¬ 
bation equations and the publicly available Boltzmann code CM- 
BEASY (Doran 2005). Initial conditions are produced using a uni¬ 
form Cartesian grid and a first-order ZeTDovich approximation us¬ 
ing a modified version of the publicly available N-GENIC code. The 
modified code uses the growth factors and expansion history calcu¬ 
lated by CMBEASY to correctly calculate the particle displacements 
in the non-standard cosmologies. 

We produce several well resolved dark matter haloes using 
hydrodynamical zoom simulations in the cosmologies mentioned 
previously. As the N-Body code used here is described in detail in 
Carlesi et al. (2014a), we briefly summarise the key points here. 
DARK-GADGET is a modified version of P-GADGET-2. The key 
modifications are the inclusion of a separate gravity tree to account 
for the additional long range forces arising from the scalar field, 
and an evolving dark matter N-body particle mass which models 
the decay of the dark matter density. The code requires the full evo¬ 
lution of the scalar field, 0, the mass of the dark matter N-body 
particle, and the expansion history. 

The zoom simulations used a parent simulation of Lbox = 
50 /i“^Mpc containing 2 x 256^ particles (DM and gas parti¬ 
cles) with the following cosmological parameters. Candidate ob¬ 
jects were identified using VELOCiraptor (Elahi et al. 2011). For 
each selected object, all particles within a radius of ~ 3i?max of 
the halo at z = 0 are identified in the high resolution IC that 
was down sampled to produce the low resolution parent simula¬ 
tion. The mass resolution in the zoom region is 8 x 10® h~^MQ 
and 1.5 X 10® h~^MQ for dark matter® and gas particles respec¬ 
tively, or an effective resolution of 1024®. All the simulations are 
started at z = 65 with the same phases in the density perturbations 
and use a gravitational softening length of 1 /40 of the interparticle 
spacing. We study our haloes across cosmic time, focusing specif¬ 
ically on z = 0 and z = 0.5, where the progenitors of the clusters 
are large galaxy groups. The cosmological parameters relative to 
the fiducial ACDM model at z = 0.5, where the quantities are not 
normalised, are listed in Table 1 . 


^ This is the dark matter particle mass at 2 ; = 0 for the coupled dark sector 
simulations. 
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Figure 1. Smoothed projected density images of the high resolution region at 2 : = 0 centred the two groups used in this study. Q (top two rows) & MM 
(bottom two rows). From left to right: ACDM, ^CDM, 0(/3 = 0.05)CDM. Top row in each subpanel shows CDM distribution and bottom row shows gas 
distribution. Densities are plotted on a logarithmic scale with bright regions being dense. The scale is arbitrary but kept fixed for a given particle type and halo. 


We identify all bound haloes and their substructures using VE- 
LOClraptor within the high resolution region. This code identifies 
field haloes using a 3D Friends-of-Friends (FOF) algorithm with 
a linking length of 0.2 times the inter-particle spacing^. Substruc¬ 
tures are then found by identifying particles that belong to a local 
velocity distribution that differs from average background and then 
linking this outlier population using a phase-space FOF approach 
(for details see Elahi et al. 2011, 2013). This technique not only 
identifies bound subhaloes, it will also identify the tidal features 
associated with the subhalo and even completely unbound tidal de¬ 
bris. 


3 HALOES IN COUPLED COSMOLOGIES 

Here we focus on two objects: one with a quiescent history 
(Q); and one which is undergoing and has undergone several major 

^ We also apply a 6DFOF to each candidate FOF halo using the velocity 
dispersion of the candidate object to clean the halo catalogue of objects spu¬ 
riously linked by artificial particle bridges, useful for disentangling mergers 
^ Weusethe Apbg = AOm,o (3if^(a)/87rG) ACDM definition, that is 
here we are not explicitly taking into account the differences in the expan¬ 
sion history when determining the overdensity that constitutes a virialised 
system in the non-standard cosmologies. 


mergers (MM). These two objects are chosen to explore the effects 
of coupled dark sector physics in two different dynamical regimes. 
We show the dark matter & gas distributions at 2 ; = 0.0 & z = 0.5 
in Fig. 1 and list their bulk properties in Table 2. 

Focusing on Fig. 1, we see that the quiescent group looks 
remarkably similar at all times in spite of the fact that (/)(/3 = 
0.05)CDM-Q halo has a maximum circular velocity that 20% 
higher and is ^ 80% more massive at z = 0 that its counterparts in 
uncoupled simulations. Considering the how self-similar dark mat¬ 
ter haloes appear, this fact in itself is not very surprising. 

MM on the other hand appears to be at different stages of a 
major merger. In = 0.05)CDM, at 2 ; = 0, the halo is relax¬ 
ing from the major merger event that occurred at z = 0.5. Three 
“cores” are visible, remnants of past mergers. The ongoing merger 
at z = 0.5 is the reason for the large differences in the total bound 
FOF mass of the halo relative to its virial mass seen in Table 2. The 
MM halo in uncoupled cosmology also has several cores at z = 0 
and is near a similar mass halo which which it will merge in the 
near future. On the other hand, ACDM-MM has yet to undergo its 
first major merger, as evidence by its lack of merger remnants. 

Clearly for our choice of cosmological parameters, haloes at 
these mass scales are more massive in strongly coupled cosmolo¬ 
gies but can have similar masses in cosmologies with the same dark 
matter density and physics but different expansion histories and 
power spectra. Hence, the increase in mass is a result of higher mat¬ 
ter densities (particularly at early times), a higher effective gravi- 
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Table 2. Bulk Properties of Q & MM: total bound mass Mtot \ total number of particles Np\ mass within the so-called virial radius = 47rApbg^A/^’ 
where A = 200; baryon fraction maximum circular velocity Vmax; radius of maximum circular velocity number of subhaloes within the virial 

radius A^subs^ mass fraction in subhaloes /M,subs^ total number of substructures (i.e., both intact and tidally disrupted subhaloes) within the virial radius Nq; 
mass fraction in substructures /m,s^ mass fraction in unbound tidal debris A/M,s,tid- 



Cosmology 

Aftot 

[1013 Mq] 

-^p,DM 

10 ® 

Ma® 

[1013 Mg] 

fh 

Arnax 

km/s 

kpc 

O [> 


Asubs /M,subs 

As 

/m,s 

A/M,S,tid 


ACDM 

8.25 

5.41 

6.44 

0.132 

687 

177 

865 

4.89 

1036 

0.069 

1050 0.098 

0.029 

Qiz = 0) 

(f)CDM 

(j){l3 = 0.05)CDM 

9.39 

16.00 

6.18 

10.73 

7.97 

11.91 

0.125 

0.119 

743 

826 

265 

314 

930 

1063 

3.51 

3.38 

1454 

2544 

0.097 

0.122 

1485 

2599 

0.135 

0.198 

0.038 

0.078 


ACDM 

4.29 

2.85 

3.40 

0.139 

629 

118 

470 

3.96 

633 

0.047 

658 

0.079 

0.032 

Q(^ = 0.5) 

0CDM 

(/)(/3 = 0.05)CDM 

4.20 

7.52 

2.76 

4.97 

3.80 

6.04 

0.126 

0.130 

673 

778 

193 

248 

486 

567 

2.51 

2.28 

529 

960 

0.071 

0.080 

542 

991 

0.102 

0.119 

0.031 

0.039 


ACDM 

3.11 

2.09 

2.94 

0.126 

518 

207 

667 

3.21 

445 

0.033 

467 

0.094 

0.059 

MM(^ = 0) 

(PCDM 

(j){l3 = 0.05)CDM 

7.04 

13.85 

4.63 

9.43 

6.39 

11.70 

0.127 

0.117 

668 

708 

393 

774 

864 

1057 

2.19 

1.36 

1434 

2293 

0.128 

0.077 

1450 

2336 

0.244 

0.113 

0.116 

0.036 


ACDM 

2.42 

1.62 

1.94 

0.115 

501 

211 

381 

1.80 

335 

0.118 

356 

0.177 

0.059 

MM(z = 0.5) 

c/)CDM 

(/)(/3 = 0.05)CDM 

2.40 

6.07 

1.61 

4.07 

2.20 

3.10 

0.123 

0.117 

553 

621 

170 

165 

398 

446 

2.33 

2.70 

315 

385 

0.068 

0.061 

335 

415 

0.142 

0.102 

0.074 

0.039 


tational force experienced by dark matter and higher growth func¬ 
tions (see Table 1). This difference is in spite of the fact that these 
parameters change by a few percent at most, and asiz = 0) is 
the same in every cosmology. However, larger haloes is not an un¬ 
ambiguous signal of a coupled dark sector. Despite the differences 
in mass at this particular scale and time, the freedom allowed in 
the (Qrn, cTs) plane means that the current differences in the mass 
of these objects is not necessarily a useful signature of coupled 
dark matter. For both Q & MM, moving from one cosmology to 
the next appears to look like viewing the same object at different 
times. Therefore, it is quite possible that a different viable choice 
of (Qrn, cTs) would remove these differences. 

None of the properties listed in Table 2 show an unambiguous 
trend with coupling that is not affected by the merger history. For 
example, the concentration parameter was found by Carlesi et al. 
(2014a) to be lower in coupled cosmology. This trend is followed 
by Q but not by MM at early times. Moreover, even the uncoupled 
cosmology has a lower concentration for these haloes. The frac¬ 
tion of the host’s mass bound up in subhaloes also varies between 
cosmologies as does the fraction in tidal debris. 

Even sampling the growth of structure across some fraction of 
cosmic time need not be necessarily effective for the same reason 
that as the observed mass difference. The weak gravitational lens- 
ing study by Pace et al. (2015) & Giocoli et al. (2015) found dif¬ 
ferences in the convergence power spectrum between ACDM and 
coupled cosmologies, however, much of the difference arises to dif¬ 
ferent as normalisations. We therefore study the internal properties 
of these objects in search of unambiguous probes of coupled dark 
sector physics. 


3.1 Profiles 

The first question is whether the density profiles differ. Carlesi et al. 
(2014b) stacked clusters to measure the average radial density pro¬ 
file at z = 0 and found negligible differences save in the poorly 
resolved, inner region for the most strongly coupled cosmology. 
Based on Fig. 1, it is not obvious major differences will exist for 
the quiescent Q group but there should be differences in the merg¬ 
ing MM group, particularly for 0(/3 = 0.05)CDM arising from 


the dynamical state. The Q quiescent halo profiles shown in Fig. 2 
are almost indistinguishable. The residuals do show that the cen¬ 
tral density of the strongly coupled cosmology is lower than the 
ACDM one but key is that the residuals are flat. The shape is not 
systematically different, all are reasonably characterised by NFW 
(or Einasto) profile (e.g Navarro et al. 1997, 2004). 

Eor MM, shown in Eig. 3, the effects of the merger are clearly 
visible. The deviation away from an NEW in the central regions of 
0(/3 = 0.05)CDM-MM at z = 0 are due to the presence of mul¬ 
tiple cores (see Eig. 1, where the cores for MM are very evident)®. 
A deviation is also seen in the uncoupled model. The profile of 
0(/3 = 0.05)CDM-MM ai z — 0.5 also extends to much larger 
radii, a result of the major merger taking place, with the smaller 
group residing just past R/\. Of greater importance is the similarity 
in the shape of the profiles at z = 0.5, when the main halo is rela¬ 
tively undisturbed in all three cosmologies. Again, the cosmologies 
have produced haloes which are remarkably similar. Even the shal¬ 
lower interior slopes in coupled cosmologies noted for Q and in 
previous studies of relaxed haloes is not seen here as it has been 
affected by the dynamical state of the system. 

Given that the effective gravitational field experienced by dark 
matter differs from that of baryons in coupled cosmologies, one 
might expect a difference in the distribution of baryons, which is 
shown in Eig. 4. Unfortunately, there is a great deal of disparity 
between the cosmologies and little indication that on an individ¬ 
ual object basis, the baryons can be used to discriminate between 
uncoupled and coupled dark matter physics, even in simplistic adia¬ 
batic simulations. By looking at various times and dynamical state, 
we can clearly see that differences in accretion history are of greater 
importance than a coupled dark sector. 


® Note the (/)(^ = 0.05)CDM-MM profile in Fig. 3 is truncated at small 
radii as a result of 1) the CM, which is calculated by an iteratively using ever 
smaller radial apertures, lies between two cores and 2) bins must contain 1% 
of the total number of particles in the halo. The CM lies off a density peak 
and as a result the first bin extend to larger radii than in the other MM halos. 
We do not attempt to address this minor issue in this figure as it emphasises 
the unrelaxed state of the halo. 
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Figure 2. Radial density profiles along with residuals (where we have re¬ 
moved the average offset to emphasise changes in slope) of Q at 2 ; = 0 (top 
two panels) and at 2 ; = 0.5 (bottom two panels). We show a thick dark gray 
region showing the best fit NFW prohle to the ACDM model and a dashed 
vertical line at the htted scale radius. 


4 SUBHALOES: COUPLED COSMOLOGIES IN HIGH 
DENSITY ENVIRONMENTS 

The bulk host haloes properties show few unambiguous differences 
between uncoupled and coupled cosmologies. Here we examine if 
subhaloes, which will host galaxies of similar mass to the Milky 
Way, show major differences. We first focus on the subhalo popu¬ 
lation residing within the high density environment inside the virial 
radius of the host haloes, Q & MM. 


4.1 Mass distribution 

The simplest comparison to make is with the subhalo mass distri¬ 
bution, shown in Fig. 5-6. Note that we do not classify the main 
merging halo in the 0(^5 = 0.05)CDM-MM simulation at 2 : = 0.5 
as a “subhalo” and so it is not present in this figure. Remarkably, 
the shape of the subhalo mass function is unaffected by the dif¬ 
ferences in the initial power spectra, expansion histories and, at 
first glance, even the dynamical state of the host halo! This re¬ 
sult is seen in the lower subpanels of each figure. The residual 
(5ni/nACDM = ni/nACDM - (ni/nACDM), has the average am¬ 
plitude difference between models, (ni/nACDM), removed to em¬ 
phasise differences in the slope. The residuals at masses above the 



Figure 3. Same as Fig. 2 but for MM. 

resolution limit, although noisy, are relatively flat or show little sys¬ 
tematic tilt, at least fox z = 0. At z = 0.5, the coupled cosmol¬ 
ogy is tilted relative to the uncoupled cosmologies. However, this 
tilt is not systematic, in one halo, the mass function is steeper, in 
the other it is flatter. Moreover, this difference is it present at even 
higher redshifts (not shown here for brevity). Clearly the shape of 
the subhalo mass distribution, which would be an ideal signature, 
is primarily governed by nonlinear dynamics (mass loss, dynami¬ 
cal friction, etc.), rather than cosmologically dependent quantities 
such as the accretion rate or the small difference in the gravitational 
force experienced by dark matter. 

The top subpanels of this figure also indicate that the differ¬ 
ences in the number of and mass fractions within subhaloes seen in 
Table 2 are not informative as the differences in the amplitude of 
the mass distribution curves are negligible and not correlated with 
the underlying cosmological model. Coupled cosmologies, at least 
for /3 ^ 0.05 are not intrinsically richer or poorer than their uncou¬ 
pled counterparts. 

4.2 Concentration 

We next examine the concentration of subhaloes in Fig. 7, specif¬ 
ically the mass dependence of ratio between the virial radius rela¬ 
tive to the maximum circular velocity radius, (we only show 

z = 0 for brevity as z = 0.5 is similar). There does appear to be a 
trend for coupled cosmologies to have less concentrated subhaloes. 


© 0000 RAS, MNRAS 000, 000-000 




















Coupled Dark Sector (Sub)Haloes 7 



Figure 4. Radial baryon fractions of Q (top two panels, with upper for 2 ; = 
0, lower for 2 : = 0.5) and MM (bottom two). Line styles and colours are 
the same as in Fig. 2. 

particularly at smaller masses, in agreement with previous studies 
(Baldi et al. 2010; Carlesi et al. 2014a). The variation in at 
a given mass is large, possibly masking this difference. Moreover, 
although some of the reduced concentration in the coupled simula¬ 
tion is due to the extra fictional force experienced by dark matter, 
some of it will also be driven by different expansion histories. The 
concentration of a halo depends in part on the formation time of 
the object (e.g. Bullock et al. 2001; Ludlow et al. 2014) and is also 
affected by the dynamical state of a (sub)halo. 

If we turn our attention more directly observable quantities, 
that is Fmax — ^Vmax Fig- the differences are less evident. 
Subhaloes are more extended at a given l/max in ACDM relative to 
(/)CDM. Although including coupling has moved the subhalo dis¬ 
tribution closer to the ACDM one, the fact that the quintessence 
model can differ from the ACDM model suggests that the concen¬ 
tration of subhaloes is not an ideal indicator of a coupled dark sec¬ 
tor. 

4.3 Angular momentum 

The angular momentum of field haloes in non-standard cosmolo¬ 
gies shows indications that there maybe subtle differences in the 
peak of the spin parameter distribution (e.g. Hellwing et al. 2013; 



Figure 5. The 2 ; = 0 cumulative mass function along (large panels) with the 
residuals in the cumulative number (smaller panels) for Q (top two panels) 
and MM (bottom two panels). In the smaller panels showing the ratio, thin 
lines correspond to where rii & tiacdm ^ 10, i.e., where the statistics are 
poor. Line styles and colours are the same as in Fig. 2 

Carlesi et al. 2014b). Here we examine the spin parameter as de¬ 
fined by Bullock et al. (2001), 

A = -;=^ -, (4) 

V2MVR 

where L is the angular momentum, M is the total mass, and R 
is the radius of the subhalo, and = GMjR. Normally, these 
quantities are measured at the so-called virial radius, i^A, however, 
subhaloes can be truncated at overdensities higher that A. For those 
subhaloes, we use the total mass, angular momentum and size of the 
object to calculate A. 

Figure 9 shows the total distribution across all mass scales. 
Note that we have calculated spins using the self-bound portions of 
a substructures, that is we do not include the loosely unbound tidal 
tails associated with subhaloes and ignore tidal debris, which would 
produce a longer high spin tail. First we note that the subhalo spin 
distribution here is lower than that found in Onions et al. (2013), 
who studied subhaloes in a galaxy mass host. The spin distribution 
at higher redshifts and lower host masses is in better agreement 
with Onions et al. (2013). 

The major feature of this plot is the absence of systematic dif¬ 
ference in the distributions between uncoupled and coupled cos- 
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Figure 6. Same as Fig. 5 but for MM. 


mologies. There is a hint in a shift towards higher spins in this cou¬ 
pled cosmology, there being a small deficit in probability relative to 
the uncoupled cosmologies just below the peak of the distribution 
and a possible excess above. However, the differences are not in¬ 
credibly significant. A simple two-sided KolmogorovSmirnov (KS) 
test comparing the distributions returns p-values of > 0.5, regard¬ 
less of which two cosmologies are compared, that is the empirical 
distributions arise from the same parent distribution. This is also 
true at higher redshifts. Moreover, the difference is more notice¬ 
able in Q than in MM, showing the dynamical state of the host 
does infiuence the spins of subhaloes residing in it. 

Even separating the distribution into masses, as we do in 
Fig. 10, does not necessarily show an unambiguous signature of 
a coupled dark sector. Small coupled cosmology subhaloes of have 
marginally higher spins. However, even the uncoupled cosmology 
differs from ACDMand for the well sampled mass scales, the mass 
dependence of A shows little difference between cosmologies with 
the two-sided KS-test giving p-values of > 0.4 regardless of mass 
bin used. Moreover, the spin parameter is infiuenced by galaxy for¬ 
mation physics, which we have not included in our simulations. 
The effect of forming stars is likely to wash out the small differ¬ 
ences noted here. 



Figure 7. The mass-concentration relation at 2 : = 0. We first bin the 
data in mass and calculate the median and the 16% & 84% quantiles of 
cvmax i^ t)in. These are plotted in the large solid points with er¬ 
ror bars. We also plot the outlying points as small points. Line styles and 
colours are the same as in Fig. 2, markers are: ACDM circle, (f)CDM square, 
0(/3 = 0.05)CDM diamond. Outlying point colours are: ACDM grey, 
(/)CDM cyan, (f){^ = 0.05)CDM orange. The x position of the bins are 
offset by an small arbitrary amount between each cosmology for clarity. As 
quantiles are misleading for a small sample, we do not plot error bars for 
bins containing fewer than 5 points, instead, we plot all the points in the bin 
and the median of these points. 



Figure 8. The Vniax — distribution similar to Fig. 7. 

4.4 Alignment 

Observations suggests that subhaloes/satellite galaxies are 
anisotropically distributed within there host halo. ACDM sim¬ 
ulations naturally give rise to anisotropic subhalo populations, 
however, this cosmology struggles to explain the strength of the 
alignments observed. Here we examine if a coupled dark sector 
produces a more anisotropically distributed subhalo population via 
the angle between a subhalo’s position, x, and major, semi-major. 
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Figure 9. Subhalo spin distribution at 2 ; = 0. For reference we also show 
the subhalo spin distribution from Onions et al (2013) for subhaloes from 
a dark matter only simulation of a 10^^ Mq halo by a solid gray line. 
Onions et al (2013) found that subhaloes peaks at a lower spin value than 
the halo distribution. Probabilities and errors are estimate using bootstrap 
resampling. Line styles and colours are the same as in Fig. 2. 



Figure 10. Spin as a function of mass. Format similar to Fig. 7. 


and minor axes of the host halo defined by the eigenvectors Ci of 
the reduced inertia tensor (Dubinski & Carlberg 1991; Allgood 
et al. 2006), 




(5) 


Here the sum is over particles in the halo, (r^)^ = (x'^Y + 
{y'n /^)^ + (^nis the ellipsoidal distance between the subhalo’s 



Od) (k2 04 06 08 LO 


cos Oi 

Figure 11. Subhalo alignment with the axes of the host halo for Q at 2 ; = 0. 
We show the probability for a isotropically distributed distribution by a gray 
line. Probabilities and errors are estimated using bootstrap resampling. Line 
styles and colours are the same as in Fig. 2. 


centre of mass and the nth particle, primed coordinates are in the 
eigenvector frame of the reduced inertia tensor, and g & s are the 
semi-major and minor axis ratios respectively. 

We show in Fig. 11 the distribution of cos Oi = x • for 
Q at 2 : = 0. This figure indicates subhaloes are typically aligned 
along the major axis of the host halo. However, haloes in coupled 
cosmologies do not have subhalo populations that are any more or 
less anisotropic than those in uncoupled cosmologies. Despite the 
vary different dynamical state of MM between cosmologies and 
relative to Q, the overall picture is similar. The subtle difference is 
that subhaloes are not as strongly aligned/anti-aligned to the ma- 
jor/minor axis of the host halo, a result of the violent phase-mixing 
in major mergers and randomisation of orbits. This type of analy¬ 
sis in itself does not indicate whether or not coupled cosmologies 
produce planar arrangement of satellites. 

To see if planes are present we examine the distributions of 
distances to the three planes with normal vectors defined by the 
inertia tensor’s eigenvectors. This distribution is compared to an 
isotropically distributed subhalo population with the same over¬ 
all radial distribution present in the halo. We show for Q the ratio 
relative to the isotropic expectation normalised by the uncertainty, 
Sp/a = (pobs/piso - l)/cr(pobs/piso), as function of distance in 

Fig. 12. 

A planar structure would show up as an excess in the proba¬ 
bility of finding subhaloes at small distances from a given plane. 
Focusing on the “major” plane, that which is perpendicular to the 
major axis that has a normal vector along the major axis, we see that 
the population is generally consistent with an isotropic distribution 
(within 2a). There is a hint of a deficit of objects at small projected 
distances that is matched by the increase in the number of subhaloes 
with small projected distances orthogonal to this plane. In general, 
< 2a deviations occur at all distances from the planes examined 
here, arising from the overall anisotropic distribution of subhaloes. 
The exception is the “minor” plane, which shows > 3a deviation 
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Figure 12. Relative probability distribution of distances to three planes nor¬ 
malised by the variance in the observed PDF for Q at 2 : = 0. The expecta¬ 
tion value for an isotropic distribution is shown by a gray line. Probabilities 
and errors are estimated using bootstrap resampling. Line styles and colours 
are the same as in Fig. 2. 


away from an isotropic distribution at small projected distances for 
all three cosmologies. 

Focusing on ACDM-Q, there are 198 subhaloes (20% of the 
subhaloes composed of ^ 100 particles) within a projected dis¬ 
tance of 40 kpc from the “minor” plane spanning a range of radii, 
from 280 kpc to 1600 kpc, that is from to past 2Ra. 

However, this planar structure is purely spatial, it is not a dynami¬ 
cally stable structure. The angular momentum of the subhaloes not 
strongly aligned with the normal vector defining the plane (45°). 
Subhalo’s are as likely to have tangentially dominated velocities as 
radial ones and be rotating, counter-rotating or have an orbit princi¬ 
pally out of the plane. The velocity dispersions perpendicular to the 
plane for even the coldest 40% is ~ 200 km/s, meaning this plane 
would disperse within 200 Myr. As expected, ACDM haloes have 
anisotropic distributions of subhaloes but no planes. 

At first glance the inclusion of coupling looks encourag¬ 
ing. The deviations away from projected distances of an isotrop¬ 
ically distributed subhalo population are more significant, although 
it is still contains similar fractions of the subhalo population as 
ACDM(?^ 20%). This planar structure has a vertical distribution 
of ~ 70 kpc, an angular extent of ~ 180°, and spans ~ 2 Mpc 
in radius. However, its dynamics again show the transient nature 
of this plane. Its angular momentum axis is misaligned with the 
plane by 70° and the vertical velocities are offset from the plane 
by 130 km/s and will disperse within ~ 190 Myr. The main com¬ 
ponent of this plane is two groups of subhaloes, which at 2 : = 0.5 
reside in the two main filaments that connect this halo’s progenitor 
to the cosmic web. 

MM, despite being at various points in a merger in the differ¬ 
ent cosmologies has similar features: anisotropic subhalo distribu¬ 
tion and no dynamically stable planar structures. This trend is true 
for both haloes at earlier redshifts as well. It is however interesting 



Figure 13. Properties of haloes around Q verses ellipsoidal distance, r'. 
Here we show the spin, A and concentration (top and bottom respec¬ 

tively). We select haloes with masses 5 x 10^/ Mq ^ M ^ 5 x 10^° Mq, 
as these quantities depend on host halo mass. Format is similar to Fig. 7. 

to note that the statistical significance of a flatten spatial distribution 
in the “minor” plane is always higher in the coupled cosmology. 

That is not to say spatially compact planes do not exist in 
other cosmologies. We note that at z = 1 both (/)CDM & 0(/3 = 
0.05)CDM Q haloes (with M ^ 3.5 x 10^^ M©) have “minor” 
planes of satellites spanning 40 kpc with an orbital angular mo¬ 
mentum axis aligned with the plane’s normal vector, although the 
motions orthogonal to the plane are large enough (200 km/s) to dis¬ 
perse the plane within a ~ 200 Myr. Despite the high redshift, these 
host masses are closer to the mass scale of the Local Group, where 
planes have been observed. 


5 FIELD HALOES: COUPLED COSMOLOGIES IN 
LOWER DENSITY ENVIRONMENTS 

The dense environment within a group appears to show little dif¬ 
ferences in the subhalo population. However it is worthwhile ex¬ 
amining whether this is true in all environments. We therefore ex¬ 
amine the bulk properties of haloes that lie outside the group as a 
function of the distance from the centre of the host. As the haloes 
are triaxial and we are interested in exploring environments at dif¬ 
ferent densities, we use the ellipsoidal distance to the halo centre, 
r = (x^ + 2/^/^^ + where coordinates are in the eigenvec¬ 

tor frame of the reduced inertia tensor, and q & s sltq the axis ratios 
(which for Q are g ~ 0.75 & s ~ 0.47), instead of the pure radial 
distance. 

Figure 13 shows the concentration and spin of small haloes 
around Q, the halo with the quiescent merger history. This fig¬ 
ure shows haloes differ between cosmologies. Focusing on the 
top panel of Fig. 13, we find haloes are more concentrated in the 
quintessence cosmology than in ACDM. Adding coupling has re¬ 
versed the trend in resulting in increasingly less concen¬ 

trated haloes. The spin distribution to show little dependence on 
on its environment but this is due to the strict removal of unbound 
particles. If we include loosely unbound particles within the FOF 
envelope of field haloes we would clearly see that haloes closer in 
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Figure 14. Baryon fraction normalised by the cosmic baryon fraction, 
fh/iS^b /^m), for Q (top) & MM (bottom). Similar to Fig. 13. 


have higher A, a result of becoming tidally disrupted. More impor¬ 
tantly, there is little difference between cosmologies in the different 
environments. These trends apply to MM, despite its active merger 
history. 

As baryons and dark matter do not feel the same gravitational 
force, we examine the baryon fractions of haloes and the depen¬ 
dence upon environment in Fig. 14 for both Q & MM. The first no¬ 
table feature is the difference between the median /b in Q and MM 
for the same cosmology. Here the dynamical state and the density 
of the local environment has affected the median /b at a given r', 
although the overall trends remain unchanged. Small haloes have 
decreasing /b < Qb/^m with decreasing radius as the hot gas 
surrounding the cluster strips haloes with increasing efficiency the 
closer a halo is to group centre. Of greater importance is the appar¬ 
ent decrease in /b when /3 > 0. The exact median baryon fraction 
depends on the dynamical state but coupling decreases the radius at 
which there are significant drop in /b. The complication lies in the 
radius at which this occurs. Q haloes show a more gradual change 
and baryon poor subhaloes appear at much larger distances than 
in MM. Additionally, the surrounding halo population in ACDM- 
MM is strongly effected the the presence of another infalling group 
mass halo, hence the sudden drop in the median /b and number of 
subhaloes containing any baryons at all. 


6 DISCUSSION & CONCLUSION 

We have studied the internal properties of group mass haloes in 
coupled Dark Energy-Dark Matter cosmologies using adiabatic 
zoom simulations. The cosmologies used in this study consist of 
a fiducial ACDM model, a quintessence model, and a strongly cou¬ 
pled DEDM model, = 0.05)CDM. In coupled models, dark 
matter decays into dark energy, giving rise to the late time ac¬ 
celerated expansion, an additional frictional force experienced by 
DM, and expansion histories and growth functions that differ from 
ACDM. It is important to note that the best fit parameters for a 
given alternative cosmology to the CMB will differ from the best 
fit parameters of ACDM. Given this freedom, we have chosen to 


pin our simulations at z = 0 so that the density parameters and ag 
match the ACDM values based on ACDM results from Planck. 

Our approach was to examine in detail two well resolved clus¬ 
ter mass haloes with different dynamical states to identify unam¬ 
biguous signatures of a coupled dark sector. We find that, consider¬ 
ing the differences between ACDM & 0CDM, the two uncoupled 
cosmologies, haloes that form in the coupled cosmologies appear 
remarkably similar to ACDM ones, even in instances where the 
total mass differs by 80%. The density profiles of our two well re¬ 
solved haloes are well characterised by an NEW (or Einasto) like 
profile when examined at times when the haloes are relaxed. The 
baryonic fraction profiles of these two haloes vary significantly but 
not in a systematic way between cosmologies. 

The surprising result is the similarity in the subhalo popula¬ 
tions of these haloes. The form of the cumulative subhalo mass 
function is seemingly unaffected by the inclusion of extra dark sec¬ 
tor physics. We do find that both subhaloes and small haloes in the 
outskirts of the group are less concentrated on average in coupled 
cosmologies. If we look at differences in the Umax — ^v^ax plane, 
which is more directly observable, we find subhaloes with small 
circular velocities are typically more extended in coupled cosmolo¬ 
gies. Unfortunately, a halo’s concentration partially depends on 
its formation time and thus depends on the choice of cosmolog¬ 
ical parameters and when a comparison between cosmologies is 
made. Eurthermore, as the concentration of subhaloes can drasti¬ 
cally change with the inclusion of galaxy formation physics, which 
we did not include, this probe of a coupled dark sector is not ideal. 

We have also examined these haloes to see if a spatially com¬ 
pact planar arrangement of satellites is present given the observed 
distributions of satellites around the only two galaxies in which 
these types of observations are possible, our Galaxy and M31 (e.g. 
Pawlowski et al. 2012; Ibata et al. 2013; Conn et al. 2013). In all 
the cosmologies studied here, we do not find strong evidence for 
a subpopulation of satellites residing in a stable co-rotating plane. 
However, we do find for both haloes at multiple epochs, our cou¬ 
pled cosmology is more likely to have a statistically significant 
“planes”. Unfortunately, these are not dynamically stable and will 
dispersing in roughly ~ 200 Myr. These dynamically unstable pla¬ 
nar structures are a natural outcome of cosmic structure formation 
(e.g. Sawala et al. 2014; Libeskind et al. 2015): Subhaloes are prin¬ 
cipally accreted along filaments. The current tension with observa¬ 
tions lies in the how compact these planes are (thicker than observa¬ 
tions suggest) and the fact that they do not have strongly correlated 
velocities. However, the increase in the statistical significance of 
planar like arrangement of subhaloes in strongly coupled groups is 
intriguing and warrants further study. Higher filamentary accretion 
should increase the likelihood of observing planar structures. Thus 
clues for a coupled dark sector may lie in the amount of material 
accreted along filaments. Whether there is an increase the observed 
velocity correlation between satellites relative to ACDM and in bet¬ 
ter agreement with observation of Ibata et al. (2014) remains to be 
seen. A complete study of the growth of filaments and the accre¬ 
tion history of haloes is beyond the scope of this paper and will be 
pursued in the future. 

We have also looked for signatures of a coupled dark sector 
in haloes residing well outside the highly nonlinear cluster envi¬ 
ronment. Although these haloes show the same trends in their bulk 
properties as a function of cosmology, i.e., lower concentrations 
in coupled cosmologies, there is little dependence on environment 
unique to coupled dark sector physics. The only environmental sig¬ 
nature noted lies in the baryon fractions: outlying haloes of a clus¬ 
ter are more likely to have lower /b than their uncoupled coun- 
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terparts at the same ellipsoidal distance or same environment. The 
distance (or density) at which there is a significant decrease in the 
average baryon fractions occurs at moderately larger radii in cou¬ 
pled cosmologies. However, this signature is subtle and is not no¬ 
ticeable when examining baryon fractions as a function of radial 
(or projected distance) distance to the group centre. Moreover, we 
have neglected all the complexity associated with galaxy formation 
physics and the exact treatment of gas and subgrid physics signif¬ 
icantly effects the ability of similar mass haloes to retain gas in a 
given environment (e.g. Tasker et al. 2008; Scannapieco et al. 2012; 
Sembolini et al. 2015). 

Thus, the internal properties of haloes are seemingly indif¬ 
ferent to a coupled dark sector, presenting an interesting possibil¬ 
ity: that moderately coupled dark sector models are viable and this 
physics could just be hidden from view. 
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